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Abstract 


Various  coirputer  models  developed  to  aid  in  the  design  of  magnetic  induction  launchers  or 
coilguns  are  described.  The  first  nwdels  the  launch  process  based  on  launcher  circuit  parameters, 
in  particular,  the  coil  inductance  as  a  function  of  projectile  position  L(x).  The  second  set  of 
models  calculates  the  overall  inductance  of  the  launcher  coil  both  with  and  without  the  plate 
present.  The  difference  between  the  values  AL  =  L(«>)  -  L(0)  has  proven  to  be  a  key  parameter 
in  determining  the  relative  performance  of  different  launcher  designs.  The  third  set  of  models 
estimate  the  fields  and  forces  generated  by  the  launcher,  under  the  assumption  of  filamentary 
conductor  structures  and  static  conditions.  The  n^thods  used  to  avoid  the  complications 
associated  with  the  dynamic  nature  of  the  actual  launch  process  are  described.  These  last  models 
yield  an  inductance  gradient  function  dL/dx  that  is  crucial  for  the  assessment  of  launcher 
performance.  The  usefulness  of  these  models  has  been  demonstrated  in  a  project  to  design  a 
coUgun  that  launches  metal  plates  against  incoming  kinetic  energy  (KE)  penetrators. 
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1.  Introduction 


The  work  reported  here  is  an  outgrowth  of  a  program  to  develop  a  system  that  can  aid  in  the 
protection  of  vehicles  from  long-rod  kinetic  energy  (KE)  projectiles.  In  the  past,  there  have  been 
investigations  into  the  possibility  of  defeating  these  projectiles  using  explosively  launched  plates 
(Frey,  Melani,  and  Stegall  1988;  Hackbarth  et  al.  1992, 1993);  however,  because  of  difficulties 
associated  with  the  use  of  explosives,  alternative  means  of  launching  these  metal  plates  are  being 
considered.  One  method  utilizes  electromagnetic  (EM)  forces  to  generate  the  large  accelerations 
needed  for  this  purpose. 

The  "reconnection  gun,"  developed  by  Cowan  et  al.  (1987, 1991)  is  one  example  of  an  EM 
plate  launcher.  This  is  a  multistage  magnetic  induction  launcher,  or  coilgun,  that  has 
successfully  launched  a  150-g  plate  at  a  velocity  of  1  km/s.  The  size  and  weight  of  a  multistage 
device  precludes  its  use  on  a  vehicle;  however,  as  a  result,  only  single-stage  devices  have  been 
considered  in  this  program  (Hummer  and  Hollandsworth  1991, 1992).  One  such  design  has  been 
used  to  test  the  effectiveness  of  metal  plates  in  defeating  subscale  metal  rods  (Hummer, 
Hollandsworth,  and  Beming  1995;  Hummer  1996).  Achieving  launch  velocities  of  several 
hundred  meters  per  second  with  a  single-stage  device  is  a  challenge  as  it  requires  internal  stresses 
that  peak  at  the  several-hundred-ton  level.  This  report  contains  descriptions  of  some  of  the 
modeling  and  simulation  tools  that  have  been  developed  to  aid  in  the  design  of  such  launchers. 
These  models  have  been  particularly  useful  in  identifying  the  key  issues  governing  the 
performance  of  induction  launchers.  The  efficacy  of  these  tools  is  demonstrated  by  comparison 
with  experiments. 

2.  Modeling  EM  Induction  Plate  Launchers 

A  magnetic  induction  launcher  is  a  simple  device  that  consists  of  an  electromagnet  coil  (or 
coils)  into  which  a  conducting  object  (the  projectile)  is  placed.  The  magnetic  forces  that  launch 
the  projectile  evolve  in  the  following  manner  a  constantly  changing  electric  current  is  passed 
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through  the  magnet  coil  so  that  a  changing  magnetic  field  is  produced  in  the  bore  of  the  launcher, 
the  changing  magnetic  field  will  cause  currents  to  flow  in  the  conducting  projectile,  and  the 
interaction  between  the  magnetic  field  and  these  "induced"  currents  will  result  in  a  magnetic 
force  being  applied  to  the  projectile,  thus  launching  it  out  of  the  coil  bore.  Induction  launchers 
are  distinguished  from  other  types  of  magnetic  launchers,  such  as  railguns,  in  that  electrical 
currents  are  not  directly  applied  to  the  projectile.  The  principle  that  current  loops  spontaneously 
form  in  the  presence  of  a  changing  magnetic  field  is  known  as  Faraday's  Law  of  Induction 
(Halliday,  Resnick,  and  Walker  1993).  This  is  the  same  principle  by  which  electrical 
transformers  operate,  and  in  fact  the  launcher  can  be  thought  of  as  the  "primary"  coil  of  a 
transformer,  and  the  projectile  itself  is  the  "secondary"  coil.  It  is  useful  to  note  that  the  induced 
currents  will  set  themselves  up  in  such  a  way  as  to  cancel  out  magnetic  fields  within  the  volume 
of  the  conducting  projectile.  While  the  principle  by  which  coilguns  operate  is  simple,  there  are 
many  launcher  and  projectile  configurations  to  choose  from,  with  no  one  being  an  obvious 
"optimal"  choice  for  a  particular  application.  The  simulation  tools  described  herein  were 
developed  specifically  to  identify  the  most  effective  design  possible  for  a  launcher  of  metal 
plates. 

hi  order  to  predict  launch  coU  effectiveness,  a  suite  of  computer  programs  has  been 
developed  that,  when  taken  as  a  whole,  serve  to  simulate  the  magnetic  induction  launch  process. 
Brief  descriptions  of  three  such  programs  are: 

(1)  The  first  simulates  the  acceleration  of  the  plate  and  associated  behavior  of  the 
launcher  circuit,  using  key  electrical  parameters  as  a  starting  point  (in  particular  the 
coil  inductance  as  a  function  of  plate  position  L(x),  which  can  be  either  measured 
directly  or  modeled  via  the  following). 

(2)  The  second  calculates  the  inductance  of  the  coil  with  the  plate  either  fully  inserted 
(L(0))  or  with  the  plate  fully  removed  (L('»)),  using  the  equation  describing  the 
inductance  of  a  straight  filament  as  a  starting  point.  The  quantity  AL  =  L('»)-L(0)  has 
proven  to  be  a  key  indicator  of  performance. 
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(3)  The  third  uses  simple  magnetostatics  to  calculate  the  magnetic  field  around,  and 
subsequently  the  magnetic  forces  on,  the  launch  coil.  This  uses  the  equation  for  the 
magnetic  field  generated  by  a  straight  filament  as  a  starting  point. 

When  a  simulacrum  of  a  conducting  plate  (also  based  on  filaments)  is  included,  the  forces  on 
the  plate  can  also  be  calculated.  From  these,  the  inductance  gradient  dL{x)ldoc  can  be  inferred. 
This  feature  bypasses  the  need  to  perform  any  complex  electrodynamics  calculations. 

These  three  programs  can  be  used  to  assess  and  compare  various  hypothetical  coil  designs. 
The  reasoning  behind  each  of  these  models  will  be  described  in  more  detail  in  the  following 
sections. 


2.1  Plate  Launch  Mechanics.  Only  an  overview  of  the  system  for  modeling  the 
acceleration  of  the  plate  will  be  presented  because  it  has  been  described  m  detail  in  previous 
publications  (Hummer  and  Hollandsworth  1991, 1992;  Hummer  1996).  The  launcher's  electrical 
circuit  is  essentially  a  series  RLC  circuit,  as  described  in  almost  any  introductory  physics  text 
(Halliday,  Resnick,  and  Walker  1993).  Applying  Kirchhoff  s  Law  for  voltages  in  a  loop  to  a 
series  RLC  circuit,  the  following  differential  equation  can  be  obtained: 

QIC  -  IR  -  d(L/)/dt  =  0,  (1) 


where  Q  is  the  charge  stored  in  the  capacitor  bank  (used  to  power  the  launcher),  C  is  the 
capacitance  of  the  bank,  I  is  the  current  in  the  circuit,  R  is  the  resistance  of  the  circuit,  L  is  the 
inductance  of  the  circuit,  and  t  is  time.  This  system  deviates  fi*om  the  simplest  RLC  circuit,  in 
that  L  is  not  constant,  but  rather  a  function  of  plate  position  x,  which  itself  is  a  function  of  time. 
There  is  also  an  additional  equation  describing  the  force  on  the  plate  during  magnetic  laimch: 


F  = 


1/2  dL 
2  6x  ’ 


(2) 
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where  dL/dx  is  the  aforementioned  inductance  gradient  function,  which  embraces  all  information 
regarding  the  magnetic  coupling  between  coils  and  plate.  For  an  existing  coil,  we  assess  the 
function  Lix)  by  simply  connecting  an  inductance  meter  and  making  measurements  while  moving 
the  plate  in  a  stepwise  fashion.  For  hypothetical  coil  designs,  computer  simulations  can  be  used 
to  approximate  this  function.  In  either  case,  a  polynomial  fit  to  the  L{x)  data  is  used  in  the 
simulation  code. 

This  picture  of  the  launch  process  is  of  course  highly  simplified.  In  an  actual  system,  the 
total  circuit  resistance  will  also  depend  on  the  plate  position,  as  the  plate  itself  is  most  definitely 
a  part  of  the  circuit.  It  in  fact  amounts  to  the  secondary  coil  in  a  transformer  circuit,  which 
decouples  from  the  primary  coil  (the  launcher)  as  the  plate  is  ejected.  In  addition,  the  use  of 
large  currents  (hundreds  of  kA)  will  lead  to  joule  heating,  and  the  high-voltage  switches  used 
(e.g.,  spark  gaps  or  ignitrons)  are  likely  to  behave  somewhat  nonlinearly.  While  more  complex 
models  do  exist  (Bernard  1987),  the  simple  model  described  previously  has  proven  to  be  quite 
sufficient  for  most  piuposes. 

Figure  1  contains  examples  of  dL/dbc  functions  for  two  entirely  different  coil  designs.  The 
'Tdox  coil"  design  is  a  six-turn  solenoid  with  a  5-in  square  cross  section,  with  the  plate  inserted 
between  the  two  center  turns  (similar  to  the  one  in  Hummer  [1996]).  The  "pancake  coil"  is 
comprised  of  two  6-in  circular  spirals  (four  turns  each),  one  above  and  one  below  the  plate.  In 
both  cases,  a  4.5-in  x  4.0-in  x  0.25-in  aluminum  plate  was  used.  Figure  1  exhibits  both  their 
respective  L{x)  functions  (as  measured  with  an  inductance  meter)  and  the  associated  dL/dx 
functions.  The  overall  change  in  inductance  (AL  =  was  0.62  pH  for  the  box  coil  and 

0.64  pH  for  the  pancake  coil. 

The  fact  that  the  dL/dx  functions  peak  at  different  locations  leads  to  a  slight  difference  in 
their  respective  coil's  implementations.  This  arises  from  the  fact  that  the  force  on  the  plate  at  any 
time  is  proportional  to  both  f  and  the  value  of  dL/dx  at  the  plate's  current  position.  The  current 
trace  1(f)  resembles  an  underdamped  sinusoid,  and  most  of  the  acceleration  occurs  during  the  first 
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Figure  1.  Inductance  and  Inductance  Gradient  vs.  Plate  Position  for  Two  Different  Coil 
Designs. 

half-cycle.  A  higher  final  plate  velocity  (Vf)  can  therefore  be  achieved  if  one  can  arrange  for  the 
plate  to  be  at  or  near  the  dL/dx  peak  position  when  I(t)  is  reaching  its  first  peak.  This  can  be 
done  through  a  careful  choice  of  the  plate's  initial  position  (Xq).  The  fact  that  the  box  coil's  dUdx 
peak  occurs  at  x  =  0,  however,  precludes  this  possibility.  As  a  result,  the  optimal  Xq  value  for  the 
box  coil  is  always  Xq  =  0.  In  the  pancake  coil  case,  where  dL/dx  peaks  at  about  x  =  4  cm,  the 
optimal  Xq  value  can  range  anywhere  firom  0.5  cm  to  3.5  cm,  depending  on  the  launch  parameters 
(plate  mass  m,  initial  bank  voltage  etc.). 

In  order  to  model  the  launch  of  a  plate,  equations  (1)  and  (2)  are  solved  using  a  fourth-order 
Runga-Kutta  technique  (Burden  and  Faires  1985).  As  the  launch  genorally  occurs  on  a 
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millisecond  time  scale,  a  1-ps  time  step  is  often  used.  For  an  existing  system,  the  actual 
capacitance  (C),  parasitic  resistance  (/?),  and  parasitic  inductance  where  L  =  +  L(x),  can 

either  be  approximated  by  fitting  a  measured  I{t)  curve  to  a  damped  sinusoid  (Hummer  and 
Hollandsworth  1991, 1992;  Hummer  1996)  or  determined  more  exactly  by  fitting  to  solutions  for 
I{t)  generated  by  this  model  (see  Figure  2).  Both  sets  of  parameters  typically  yield  essentially 
equivalent  predictions  for  the  final  plate  velocity,  which  in  turn  tend  to  closely  agree  with 
measured  values.  This  is  illustrated  in  Table  1,  where  actual  and  predicted  velocities  are 
compared  for  three  shots  of  the  box  coil. 


Figure  2.  A  Typical  Fit  of  Hate  Launch  Model  Resulte  to  Measured  Current  Data. 

The  plate  launch  model  is  useful  for  predicting  and  comparing  the  actual  effectiveness  of 
various  coil  designs.  As  a  case  in  point,  one  might  hypothesize  that  the  pancake  coil  would  be  a 
better  design  than  the  box  coil,  because  one  can  "tune"  Figure  3  shows  that  this  is  not  die 
case,  in  that  simulations  predict  almost  identical  v/s  for  both  coils  (for  any  given  initial  voltage 
Vq)'  The  parameters  used  in  these  simulations  were  C  •=  1,600  pF,  =  2.0  pH,  R  =  15.0  m  £2, 
and  m  -  200  g.  Predicted  peak  currents  (a  measure  of  peak  stress  on  the  coil)  are  also  effectively 
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Table  1.  A  Comparison  of  Predicted  vs.  Actual  Final  Velocities  for  Three  Shots  of  the  Box 
Coil.  Method  No.  1  Uses  R,  L^y^,  and  C  Values  Derived  from  a  Fit  of  I(t)  to  a 
Damped  Sinusoid.  Method  No.  2  Uses  Values  Derived  from  a  Fit  to  Functions 
Generated  by  the  Model  Described  in  the  Text 


Shot 

No. 

^0 

(kV) 

n 

Actual  Vf 
(m/s) 

Vf  pred.  (m/s) 
method  no.  1 

Vf  pred.  (m/s) 
method  no.  2 

15 

14.00 

208.14 

185 

184.5 

186.7 

16 

13.99 

206.87 

186 

185.2 

186.5 

17 

13.99 

217.33 

176 

177.8 

179.2 

Figure  3.  Theoretical  Final  Plate  Velocities  vs.  Initial  Voltage  for  Two  Different  Cofl 
Designs. 


identical.  The  nearly  identical  behavior  of  these  coils  is  due  to  their  nearly  identical  AL's.  The 
shapes  of  their  L(x)  and  dL/dx  curves  (Figure  1)  seem  to  be  relatively  unimportant  in  determining 
their  overall  effectiveness. 
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The  ability  to  simulate  the  launch  process  is  essential  if  one  wishes  to  compare  the 
performance  for  different  coil  geometries  and  launch  parameters.  Table  2  contains  the  results  of 
a  comparison  of  four  designs,  including  a  hypothetical  two-stage  device  based  on  round  pancake 
coils.  The  triggering  of  the  second  stage  was  optimized  in  order  to  achieve  the  stated  velocity 
and  energy  efficiency.  Also  included  is  a  design  based  on  squared  spiral  pancake  coils,  which 
proves  to  be  twice  as  efficient  as  the  round  pancake  coil  case,  even  though  a  plate  with  twice  the 
mass  is  used.  ITie  parameters  used  in  each  simulation  were  Vq,  =  20  kV,  C  =  1,600  pF, 

=  1.5  pH,  and  /?  =  15.0  m  Q. 


Table  2.  A  Comparison  of  Four  Different  Induction  Launcher  Designs  Based  on 
Simulations 


CoU  type 

w(g) 

y,(m/s) 

eff.  (%) 

1^0^) 

box  coil 

200 

350 

3.8 

331 

round  pancake 

200 

338 

3.6 

335 

2  stage  round  pancake 

200 

432 

5.8 

333 

square  pancake 

400 

337 

7.1 

277 

2.2  Inductance  Calculations.  As  the  parameter  AL  -  L(~)  -  L(0)  has  proven  to  be  a  good 
indicator  of  overall  launcher  performance,  it  would  be  useful  to  be  able  to  predict  these  values. 
The  calculation  of  the  inductance  of  the  launcher  coils  without  the  plate  present  L(~)  is  fairly 
straightforward.  The  coils  are  divided  up  into  convenient  segments,  the  self-inductance  of  each 
segment  is  calculated,  and  then  the  mutual  inductance  between  each  pmr  of  segments  is 
calculated.  The  total  inductance  of  the  coils  is  then  simply  the  sum  of  all  the  self-inductances 
and  mutual  inductances.  The  situation  when  the  plate  is  inside  the  coil  L(0)  proves  to  be 
somewhat  more  complex  because  the  mutual  inductance  between  each  segment  and  the  entire 
plate  must  also  be  calculated.  If  the  function  dL/dx  is  available,  it  is  much  more  convenient  to 
calculate  L(0)  using 

L(0)  =  L(<»)  -  j  (dL/dx)dx.  (3) 

0 

A  simple  method  for  calculating  dL/dx  is  given  in  a  later  section. 
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2.2.1  Calculating  Self-inductance  With  No  Plate  Present.  Detailed  descriptions  of  methods 
for  calculating  inductances  can  be  found  in  Grover  (1946)  and.,  to  a  lesser  extent,  in  Knoepfel 
(1970).  As  both  books  are  now  out  of  print,  a  synopsis  will  be  provided  here.  The  discussion 
will  be  restricted  to  coil  structures  that  consist  of  a  series  of  long,  straight  conductors  that  run 
either  parallel  or  perpendicular  to  one  another.  These  are  sufficient  to  describe,  or  at  least 
approximate,  box  coil  and  square  pancake  coil  designs.  The  latter  is  similar  to  the  round  pancake 
coils  described  earlier,  except  that  a  squared  spiral  pattern  is  milled  into  the  conductor  plate. 

Each  segment  of  the  coil  structure  is  assumed  to  be  a  long,  straight  conductor  of  some  length 
/,  with  arbitrarily  shaped  cross  section  having  area  A.  For  the  purpose  of  calculating  the 
segment's  self-inductance,  its  cross  section  is  divided  into  elemental  areas  (dA  =  dxdy),  each 
representing  a  long,  straight  filament  of  the  same  length  as  the  segment.  The  self-inductance  of 
the  segment  is  assumed  to  be  the  sum  of  the  mutual  inductances  between  all  possible  pairs  of 
these  parallel  filaments.  For  a  pair  of  parallel  filaments,  one  located  at  (x,  y)  and  the  other  at 
(x\  y’),  the  exact  mutual  inductance  is  given  by 


0.002/ 
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1 

—  + 

1  d 

1  +  —  +  — 

[d  ^ 

/2  l\ 

dxdydx  'dy 


(4) 


where  /  is  the  length  of  each  segment,  d  is  the  distance  between  the  segments,  the  distances  are  in 
centimeters,  and  the  inductance  is  in  microhenries.  If  the  distance  between  the  two  filaments  is 
much  smaller  than  the  length,  this  equation  can  be  approximated  by 


dL^  =  0.002/ 


21 

Cn(4)  -  1 

a 


dxdydx  'dy 


(5) 


This  can  be  integrated  to  find  the  total  self-inductance  of  the  segment 


=  0.002/  [«n(2/)  -  1] 


A2 


fill 


(6) 
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In  keeping  with  convention,  the  last  integral  is  related  to  a  "geometric  mean  distance"  p 
between  all  filaments  covering  the  area  A,  which  is  defined  by 


fn 


"  A  /  /  /  /  ^^”  dldcdydx^dy', 


(7) 


so  that  the  self-inductance  of  the  segment  can  be  written  as 


=  0.002/ 


Cn(2-)  -  1 

P 


(8) 


for  any  long  conductor  with  any  cross-sectional  area,  as  long  as  the  mean  geometric  distance 
across  that  area  is  known.  The  windings  of  box  coils  and  square  pancake  coils  have  simple 
rectangular  cross  sections,  and  the  integral  in  equation  (8)  can  be  expressed  in  closed  form.  The 
form  of  this  equation  is  somewhat  cumbersome,  however,  and  it  is  far  more  convenient  to  use  the 
approximation  p  =  0.223  (W  +  if),  where  W  and  H  are  the  width  and  height  of  the  rectangles, 
respectively.  This  equation  was  determined  empirically;  for  example,  one  can  calculate  that 
p  =  0.22313  (W)  for  a  thin  sheet  (H  =  0),  and  p  =  0.22353  (2W)  for  a  square  (H  =  W).  Using  the 
approximate  relation,  it  is  a  simple  matter  to  calculate  the  total  self-inductance  of  each  segment. 
The  total  self-inductance  of  the  coil  is  then  simply  the  total  of  all  of  the  segments'  self¬ 
inductances. 


2^2  Calculating  Mutual  Inductance  With  No  Plate  Present  One  can  express  the  mutual 
inductance  between  parallel  segments  of  equal  length  is  in  a  similar  fashion  as  before  (once  again 
assuming  that  distances  are  in  centimeters  and  inductances  are  in  microhenries); 


-'M 


=  ±0.002/ 


■ 

-  1 

fin 

(9) 
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where  /  is  the  length  of  the  segments,  and,  here,  is  the  geometric  mean  distance  between 
elements  of  the  area  of  the  first  segment  (AJ  and  elements  of  the  area  of  the  second  segment  (Aj). 
Expressing  die  locations  of  arbitrary  examples  of  these  elements  as  (x^,  yj  and  (x^,  yj,  p^,  is 
then  defined  by 

Iff! 

where  d^i,  is  the  distance  between  the  area  element  in  segment  "a"  and  the  area  element  in 
segment  "b."  The  integral  is  analytic  if  both  cross-sectional  areas  are  rectangular;  however,  the 
final  expression  is  once  again  rather  complicated.  In  practice,  is  well  approximated  by  the 
distance  between  the  centers  of  the  two  rectangles,  even  when  the  segments  are  quite  close 
together.  Essentially,  the  total  mutual  inductance  is  calculated  as  if  all  of  the  coil  segments  were 
replaced  by  filaments  located  at  their  centers;  the  mutual  inductances  between  all  pairs  of 
filaments  are  calculated  and  added  together.  The  only  difference  between  this  and  the  self¬ 
inductance  calculation  is  that  the  relative  direction  of  current  flow  is  now  important:  the  term  in 
equation  (9)  is  positive  if  the  currents  are  parallel,  and  negative  if  they  are  antiparallel.  Pairs  of 
segments  in  which  the  currents  flow  perpendicular  to  one  another  are  of  no  concern,  as  it  can  be 
shown  that  their  mutual  inductances  are  zero. 

As  the  segments  in  a  square  pancake  coil  are  of  unequal  lengths,  it  behooves  us  to  deal  with 
this  situation.  Unequal-length  filament  cases  can  be  broken  down  in  a  combination  of  several 
equal-length  cases.  For  example,  picture  two  parallel  filaments  with  lengths  I  and  m  (where 
l>m)  and  separation  d.  Furfliermore,  assume  that  the  longer  filament  extends  past  the  shorter  a 
distance p  at  one  end  and  a  distance  q  at  the  other  (i.e.,  l  =  m+p  +  q).  The  mutual  inductances 
between  these  filaments  can  then  be  calculated  from  four  equal-length  case  terms; 

where  the  subscripts  draiote  the  filament  length  to  be  used  in  equation  (9).  The  separation  is  in 
all  four  cases.  There  are  similar  rules  for  cases  in  which  the  segments  overlap  differently  or  not 
at  aU  (Grover  1946). 
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2.23  Results.  Two  aluminum  pancake  coil  models  were  constructed  in  order  to  test  this 
system,  each  consisted  of  two  squared  spiral  coils  with  a  ccnter-to-center  separation  of  2  in.  The 
first  coil  was  constructed  by  milling  out  a  1/4-m  wide  square  spiral  pattern  in  a  6-in  square  plate 
of  1/4-in-thick  aluminum,  leaving  behind  1/4-in  square  conductors.  There  were  five  turns  in 
each  coil.  The  second  was  constructed  similarly,  except  that  1/2-in  conductors  and  spaces  were 
formed  for  a  total  of  two  and  a  half  turns  in  each  coil.  In  both  cases,  all  inside  comers  in  the 
conductors  were  left  with  a  round  fillet  with  a  radius  equal  to  that  of  the  milling  bit  used  (a 
feature  not  included  in  the  calculations).  Several  of  the  conductor  segments  are  not  long 
compared  to  their  widths,  in  which  case  the  validity  of  equation  (8)  comes  into  question.  Even 
so,  the  agreement  was  excellent:  the  total  inductance  of  the  five-tum  coils  was  calculated  to  be 
5.92  pH  and  measured  to  be  5.91  pH.  In  the  two-and-a-half-tum  case  it  was  1.48  pH 
(calculated)  vs.  1 .54  pH  (measured). 

23.4  Calculating  Inductances  Wdh  a  Plate  Present.  The  scheme  for  estimating  the 
inductance  of  these  coils  with  a  metal  plate  inserted  is  somewhat  cumbersome,  and  only  a  brief 
description  will  be  provided.  The  difficulty  is  in  modeling  the  plate  itself,  as  it  is  a  conductor 
whose  length  and  width  are  comparable  and  in  which  eddy  currents  flow  in  a  complex  pattern. 
The  scheme  relies  on  Lenz's  Law  (Halliday,  Resnick,  and  Walker  1993),  which  implies  that  the 
eddy  currents  will  form  in  such  a  way  as  to  largely  cancel  the  field  within  the  volume  of  the  plate. 
Another  way  of  looking  at  this  is  that  the  plate  tends  to  shield  the  top  coil  from  the  bottom  (it 
would  isolate  them  completely  if  the  plate  covered  an  area  much  larger  than  the  coils 
themselves).  So,  even  though  the  plates  actually  only  cover  the  same  area  as  the  coils  (and  thus 
some  leakage  around  the  edges  is  likely),  it  is  assumed  that  there  is  no  mutual  inductance 
between  any  of  the  upper  and  lower  coil  segments,  and  these  terms  are  simply  excluded  from  the 
previous  calculation. 

Each  of  the  two  coils  does  have  mutual  inductance  with  the  plate,  however,  and  this  must  be 
accounted  for.  In  order  to  mimic  the  cancellation  of  field  at  the  plate  surface,  an  "image  coil"  is 
constructed  on  the  opposite  side  of  the  launcher's  midplane.  This  image  coil  is  identical  to  the 
real  coil  in  all  ways,  except  that  the  current  flows  throu^  it  in  the  qiposite  direction,  and,  thus. 
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symmetry  guarantees  that  the  fields  will  cancel  at  their  midplane.  The  mutual  inductance  of  the 
coil/plate  combination  is  then  assumed  to  be  the  same  as  that  for  the  coil/image  coil  combination. 
An  identical  procedure  is  then  applied  to  the  other  coil,  and  all  terms  are  added. 

The  results  for  this  calculation  are  not  quite  as  good  as  in  the  previous  case:  the  inductance 
calculated  for  the  five-turn  coils  was  3.5 1  pH,  as  compared  to  a  measured  value  of  3.58  pH.  For 
the  two-and-a-half-tum  case  it  was  0.85  pH  (calculated)  vs.  1.01  pH  (measured).  Possible 
sources  of  error  include  imperfect  isolation  of  the  upper  and  lower  coils  and  the  fact  that  the 
plate's  self-inductance  is  not  accounted  for.  This  model  has  proven  to  be  of  limited  usefulness 
because  it  carmot  easily  be  extended  to  cases  where  the  plate  offers  only  partial  shielding.  The 
following  model  offers  a  way  to  gain  this  information,  while  also  yielding  estimates  of  the 
iutracoil  stresses  and  external  magnetic  fields  and,  so,  will  be  described  in  more  detail. 

23  Quasi-static  Model.  The  launch  of  a  plate  by  magnetic  induction  is  a  highly  dynamic 
process:  the  current  in  the  coil  is  changing  in  time;  the  magnetic  field  and  eddy  currents  induced 
in  the  plate  change  along  with  it;  the  magnetic  field  attempts  to  diffuse  into  the  edges  of  all 
conductors;  the  temperature  and,  thus,  resistivity  of  the  conductors  rise,  etc.  Most  of  these  events 
can  be  simulated  by  the  application  of  Maxwell's  equations,  which  govern  all  classical  EM 
behavior.  In  practice,  however,  it  is  extremely  difficult  to  apply  the  equations  to  arbitrary  three- 
dimensional  situations.  This  was  the  motivation  behind  the  conception  of  the  following  model, 
first  described  in  (Beming  et  al.  1997),  which  conveniently  ignores  much  of  the  three- 
dimensional  nature  of  the  problem  and  almost  all  of  its  dynamic  nature. 

23.1  Magnetic  Field  and  Lorentz  Force  Calculations.  The  primary  function  this  model  is 
to  estimate  the  magnetic  field  B  created  by  the  coil  assembly.  Considering  the  case  where  no 
plate  is  present,  we  begin  by  assuming  that  the  magnitude  of  the  field  (at  any  point)  at  a  particular 
time  is  proportional  to  die  current  in  die  coil  at  that  time.  As  such,  we  completely  ignore  delayed 
responses  of  the  system.  It  also  implies  that  no  ferromagnetic  material  is  present  in  the  vicinity. 
This  assumption  allows  us  to  calculate  the  field  using  magnetostatics.  Furthermore,  we  exploit 
the  fact  that  the  B-field  outside  of  a  length  of  conductor  of  finite  cross  section  tends  to  resemble 
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the  field  due  to  an  infinitesimally  thin  filament  of  the  same  length,  which  is  simple  to  describe. 
For  a  finite  length  filament  carrying  current  I  and  oriented  parallel  to  the  x-axis,  having  endpoint 
coordinates  (xoi,  yo,  Zq)  >’o»  di®  magnetic  field  components  at  an  arbitrary  point  in 

space  (x,  y,  z)  are,  in  MKS  units. 


5, 


r  \ 

,  An  j 


(z  -  Zq) 


(Xo2  -  x)r^-  (Xqj  -  x)rj  -  (rj)' 
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t  T  \ 

47T  ) 


(y  -  >o) 


(Xo2  -  X)r2  -  ir^)^  (Xoi  -  x)rj  -  (r^y 


where 


=  )j(x  -  Xoj)2  +  (y  -  y^y  +  (z  -  Zq)^  , 


and 

^2  =  ^(x  -  Xo2)^  +  (y  - 


Alternate  orientations  would  require  the  appropriate  coordinate  transformation.  One  practical 
concern  involving  the  equation  (12)  is  that  it  yields  “<»”  at  all  points  along  the  segment  itself. 

This  is  due  to  the  unphysical  assumption  of  zero  thickness  and  can  be  alleviated  in  a  simulation 
program  by  forcing  B  to  zero  within  some  arbitrarily  tiny  distance  around  the  segment  itself. 

The  program  “BFELDIO"  has  the  ability  to  simulate  coil  structures  with  combinations  of 
filamentary  segments,  provided  each  is  oriented  in  either  the  x,  y,  or  z  directions.  Currently,  it 
can  automatically  construct  pancake  coils  consisting  of  concentric  rectangles  or  rectangular 
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spirals,  given  the  required  dimensions.  It  can  also  approximate  the  stmcture  of  a  box  coil  with 
coaxial  rectangular  loops.  It  then  calculates  magnetic  field  strength  and  direction  and  can  display 
these  graphically.  It  also  estimates  the  magnetic  forces  felt  by  various  sections  of  the  launcher  by 
dividing  each  filament  into  a  number  of  segments  and  then  calculating  the  Lorentz  force  F  felt  by 
that  segment  The  Lorentz  force  is  given  by  (Halliday,  Resnick,  and  Walker  1993): 

F^ILxB,  (13) 

where  L  is  the  length  of  the  segment,  and  B  is  the  field  at  the  segment's  center.  These  results  can 
also  be  displayed  graphically. 

232  Plate  Models  and  Inductance  Gradient  Calculations.  The  aforementioned  system 
treats  the  coil  assembly  as  if  it  were  a  simple  DC  electromagnet  This  makes  it  difficult  to  model 
the  plate  because  all  of  the  currents  present  in  the  plate  during  launch  are  a  direct  result  of  the 
changes  in  the  magnetic  field.  Lenz's  Law  implies  that  the  eddy  currents  induced  in  the  plate  will 
take  forms  that  attempt  to  cancel  the  field  within  the  plate  itself  and  will  constantly  adjust  to 
follow  the  external  field.  How  well  these  currents  succeed  in  expelling  the  field  depends  on  the 
conductivity  of  the  plate;  a  perfect  conductor  would  have  zero  field  throughout  its  volume,  and  a 
less  than  perfect  conductor  would  allow  some  field  in  along  its  outer  surfaces.  The  characteristic 
distance  of  penetration  in  referred  to  as  the  "skin  depth,"  and  is  typically  less  than  2  mm  for  the 
metals  and  firequencies  used  in  this  study  (Knoepfel  1970).  The  actual  eddy  currents  that  perform 
this  feat  clearly  must  possess  a  highly  detailed,  three-dimensional  structure,  and  determining  that 
structure  would  require  a  difficult  three-dimensional  electrodynamics  calculation. 

The  following  model,  which  uses  a  finite  number  of  filaments  to  model  the  plate,  must 
therefore  forego  any  attempt  at  an  exact  simulation.  Furthermore,  the  complex  electrodynamic 
process  is  in  no  way  simulated,  but  rather,  the  effect  is.  The  scheme  is  as  follows:  (1)  a  coil 
design  and  current  are  chosen,  and  its  magnetic  field  is  simulated  previously  mentioned;  (2)  a 
basic  structure  of  filamentary  loops  is  developed  in  order  to  simulate  the  plate;  then  (3)  the 
currents  within  these  "plate"  filaments  are  adjusted  in  such  a  way  as  to  minimize  (not  cancel)  the 
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magnetic  field  within  the  plate.  This  last  step  is  a  simply  an  attempt  to  mimic  the  expulsion  of 
magnetic  field  from  the  plate.  The  coarseness  of  a  filamentary  loop  structure  precludes  an  exact 
cancellation,  so  we  are  forced  to  resort  to  some  sort  of  minimization  scheme. 

The  question  now  becomes:  what  is  the  simplest  combination  of  plate  structure  and 
minimization  scheme  that  will  yield  an  adequate  approximation  of  the  situation?  The  method  for 
determining  "adequacy"  is  quite  simple;  once  the  currents  in  them  are  set,  the  magnetic  forces  on 
the  plate  filaments  can  be  calculated  in  the  same  fashion  as  for  the  coil  filaments  and  then  added 
in  order  to  find  the  net  magnetic  force  on  the  plate.  If  the  expulsion  of  the  field  has  been 
adequately  approximated,  then  the  value  of  dL/dx  can  be  inferred  using  equation  (2),  even  though 
no  proper  inductance  calculation  has  been  performed.  Functions  derived  in  this  manner  are  then 
compared  to  measured  functions  in  order  to  evaluate  the  particular  scheme  used. 

In  order  to  illustrate  the  level  of  complexity  required  by  this  model,  a  brief  history  of  its 
evolution  follows.  The  first  version  used  a  single  rectangular  loop  covering  the  same  area  as  the 
plate,  with  its  current  adjusted  so  that  the  field  at  the  plate's  central  point  is  zero.  When  tested  on 
the  square  pancake  design  mentioned  previously,  this  yielded  dzydx  curves  that  were  poorly 
shaped  and  far  too  high.  Recognizing  the  fact  that  the  coil's  field  tends  to  have  different 
directions  inside  and  outside  the  coil  the  next  plate  model  broke  the  plate  into  two  loops:  one 
covering  the  area  of  the  plate  inside  the  coil  and  one  covering  the  area  of  the  plate  outside,  so  that 
these  regions  would  be  handled  separately.  When  the  two  plate-loop  currents  were  adjusted  to 
cancel  the  field  at  their  respective  loops’  central  points,  the  resulting  dL/dx  curves  took  on  an 
excellent  shape,  but  were  still  at  least  a  factor  of  3  too  high.  At  this  point,  a  new  minimization 
scheme  was  tried:  forcing  the  average  of  the  field's  vertical  components  inside  a  particular  loop 
to  go  to  zero.  A  finite  number  of  points  within  the  loop  were  used  in  order  to  determine  the 
average,  and,  at  first,  25  equally  spaced  points  were  used.  This  inraiediately  yielded  a  near 
perfect  dL/dx  curve  for  the  square  pancake  model,  with  the  same  excellent  shape  and 
normalization  about  10%  off.  Splitting  each  loop  into  an  “upper”  and  “lower”  one,  in  order  to 
mimic  the  three-dimensional  nature  of  the  plate,  improved  the  agreement  even  further. 
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Unfortunately,  this  remarkable  agreement  proved  to  be  partially  the  result  of  a  fortuitous 
choice  of  sampled  points,  as  the  normalization  of  the  inferred  dL/dx  curve  was  later  shown  to  be 
highly  sensitive  to  the  number  and  spacing  of  said  points,  hi  fact,  the  spacing  of  the  25  points 
can  arbitrarily  be  adjusted  so  that  the  inferred  curve  is  a  near  perfect  fit  to  the  data,  at  least  in  this 
one  case.  Using  larger  numbers  of  points  tended  to  degrade  the  fit.  The  presence  of  arbitrary 
parameters  (i.e.,  the  number  and  spacing  of  sampled  points)  led  to  the  abandonment  of  this 
scheme.  It  is  instructive,  however,  to  note  just  how  little  filamentary  structure  is  needed  in  order 
to  obtain  reasonable  semiquantitative  results. 

While  the  four-loop  model  gave  reasonable  results  using  the  minimization  of  average  field 
scheme,  it  clearly  lacked  sufficient  stmcture  to  mimic  the  actual  cancellation  of  the  complex  coil 
field.  For  this  reason,  more  structure  was  added  in  the  final  version  of  the  model;  the  plane  of 
the  plate  is  divided  into  100  equally  sized  rectangular  areas,  and  two  rectangular  loops  (an  upper 
and  lower)  are  used  in  each  region,  as  illustrated  in  Figure  4.  The  current  in  each  loop  pair  is 
adjusted  so  that  the  vertical  component  of  the  B  field  at  the  central  point  of  their  volume  is  zero. 
In  practice,  the  plate  currents  are  adjusted  so  that  their  fields  are  opposite  what  is  produced  by  the 
coil.  As  the  plate's  field  now  has  contributions  from  each  of  the  100  plate  loop  pairs,  adjusting 
the  100  plate  loop  currents  to  simultaneously  zero,  the  field  at  the  100  points  becomes  somewhat 
complicated.  However,  since  each  loop's  contribution  is  proportional  to  that  loop's  current,  and 
the  total  is  just  the  sum  of  these  contributions,  the  problem  lends  itself  to  linear  algebra.  As  this 
problem  requires  the  reduction  of  a  100  x  100  matrix,  a  highly  efficient  numerical  method  (such 
as  "LU  decomposition"  [Press  et.  al  1986])  is  recommended.  This  feature  is  included  in  the 
program  BFELDIO. 

Figure  5  illustrates  the  effectiveness  of  the  200-loop  model;  it  compares  BFIELDIO  results 
with  measured  results  for  the  6-in  x  6-in  square  pancake  described  previously.  The  coil-center  to 
coil-center  spacing  was  2  in,  and  the  plate  was  a  6-in  x  6-m  x  0.25-in  piece  of  aluminum.  The 
dL/dx  results  for  a  aude  replica  of  the  coil,  made  from  thin  copper  wire,  are  also  shown.  The 
close  resemblance  of  the  replica's  dL/dx  curve  to  that  for  the  actual  coil  is  typical,  a  fact  that 
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Figure  4.  Filamentary  Structures  Used  to  Model  a  Square  Plate  and  Square  Pancake  Coil 
System. 


Figure  5.  Measured  Inductance  Gradient  of  a  6-in  Square  Pancake  Coil  Compared  to  the 
Results  of  Simulations  Performed  by  the  Program  BFIELDIO  (for  Both  the 
Actual  6-in  Square  Plate  and  a  Hypothetical  5.5-in  Square  Hate).  Also  Included 
Are  Data  from  a  Narrow  Gauge  Wire  Replica  of  the  Cml. 


makes  replicas  a  useful,  inexpensive  method  of  testing  different  coil  designs  (although  their  L(x) 
curves  differ  by  a  constant,  due  to  the  thin  wire's  large  self-inductance).  As  seen  Figure  5,  the 
simulated  dL/dx  ciuwe  is  slightly  higher  than  and  offset  from  the  measured  curves.  This  offset 
was  not  seen  in  the  previous  four-loop  model  and  is  suspected  to  be  an  artifact  associated  with 
this  particular  plate  model.  Simulations  using  lesser  numbers  of  plate  loops  (50, 72, 98, 128,  and 
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162)  show  a  reasonable  convergence  toward  the  200-loop  answer  and  indicate  that  more  loops 
would  not  improve  the  agreement  significantly. 

Experience  indicated  that  both  the  normalization  and  the  offset  could  be  improved  if  one 
assumed  an  "effective"  plate  size  smaller  than  the  actual  size.  Figure  5  includes  results  of 
simulations  in  which  a  5.5-in  x  5.5-in  plate  size  was  assumed,  and  these  are  indeed  a  better  fit  to 
the  measured  results.  This  improvement  is  also  seen  in  Figure  6,  where  simulated  and  actual 
results  for  a  10-in  x  6-in  square  pancake  coil  and  12-in  x  6-in  plate  are  found,  assuming  an 
1 1.5-in  X  5.5-in  plate  once  again  improved  the  fit.  It  would  be  tempting  to  explain  the  efficacy  of 
"discounting"  the  plate's  edges  by  invoking  some  sort  of  skin-depth  argument;  however,  the 
0.25-m  adjustments  to  each  edge  are  much  larger  than  the  <0.08-in  skin  depth  expected  here.  It 
may  be  an  indication  that  the  situation  near  the  edges  is  not  being  properly  treated,  however,  and 
perhaps  a  finer  mesh  there  would  improve  the  fit.  Regardless,  the  unadulterated  200-loop  model, 
for  all  its  crudity,  is  still  a  remarkably  effective  device. 


X  (cm) 


Figure  6.  Measured  Inductance  Gradient  of  a  10>in  x  6-in  Rectangular  Pancake  Coil 

Model  Compared  to  the  Results  of  Simulations  (for  Botii  the  Actual  12-in  x  6-in 
Plate  and  an  11^-in  x  5^-in  Plate). 
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2.3.3  Magnetic  Field  Results.  The  magnetic  fields  generated  by  this  type  of  plate  launcher 
have  been  the  objects  of  two  prior  studies;  one  concerning  a  round  pancake  coil  (Cobum,  Le,  and 
Martin  1995)  and  the  other  dealing  with  a  square  pancake  coil  design  (Beming  et  al.  1997).  As 
the  quasi-static  model  was  utilized  in  the  latter  study,  a  brief  discussion  of  the  outcome  will  be 
included  here. 

The  launcher  in  question  consisted  of  two  6-in-wide,  4  1/2-tum  square  pancake  coils, 
machined  from  a  copper-beryllium  alloy.  The  conductors  were  0.25  in  wide  and  0.5  in  high,  with 
0.25  in  of  space  between  them.  The  coils  were  completely  encased  in  a  glass-epoxy  composite 
(GIO)  support  stmeture,  and  14-in-long  sections  of  steel  U-chaimel  were  added  to  the  top  and 
bottom  in  order  to  stiffen  the  structure.  The  aluminum  alloy  plates  were  6  in  square  x  0.25  in 
thick  and  had  a  mass  of  400  g.  The  launcher  was  fiired  many  times  from  a  small  test  stand, 
depicted  in  Figure  7,  which  contained  a  100-kJ  capacitor  bank. 


Figure  7.  The  Test  Stand  Used  in  a  Study  of  the  EM  Fldds  Generated  hy  an  Induction 
Launcher. 
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The  three  components  of  magnetic  field  were  measured  at  10  separate  points,  some  along  the 
shotline,  some  beside  the  launcher.  As  expected,  the  time  evolution  of  the  magnetic  field 
followed  that  of  the  launcher  current  Figure  8  contains  a  plot  of  peak  field  strength  vs.  radial 
distance  from  the  center  of  the  launcher  for  all  points.  The  complex  structure  seen  in  Figure  8 
would  seem  to  belie  the  simple  power-law  response  predicted  by  the  quasi-static  model  (depicted 
by  the  solid  line);  however,  inclusion  of  the  power-feed  cable  in  the  model  recreates  most  of  the 
structure.  The  power-feed  cable  was  clearly  the  dominant  source  of  field  in  some  locations. 

The  level  of  agreement  between  simulation  and  reality  seen  in  Figure  8  would  appear  to  be 
quite  reasonable,  especially  since  the  simulations  do  not  include  the  effects  of  shielding  by  the 
various  metallic  structures  in  and  around  the  launcher  (except  by  the  plate),  nor  are  all  current 
paths  included. 

23.4.  Intracoil  Stresses.  The  quasi-static  model  can  also  estimate  the  forces  on  sections  of 
the  coil  itself,  as  illustrated  in  Figure  9,  in  which  forces  acting  on  select  parts  of  a  square  pancake 
coil  are  illustrated.  Forces  are  highest  in  the  region  behind  the  plate,  as  the  field  is  concentrated 
by  the  exclusion  effect.  Forces  are  also  concentrated  at  comers,  with  total  magnitudes  on  the 
order  of  140  ksi  at  the  300-kA  level.  This  is  near  the  155-ksi  yield-strength  limit  for  the  copper- 
beryllium  alloy  used  to  constract  a  launcher  much  like  the  one  depicted.  Based  on  this 
simulation,  one  might  therefore  conclude  that  this  laimcher  could  survive  a  300-kA  pulse; 
however,  the  situation  is  much  more  complex  than  this.  In  fact,  this  launcher  failed  (via  cracks  at 
two  rear  comers)  at  a  peak  current  level  of  approximately  200  kA,  at  which  point  the  internal 
stresses  are  less  than  half  that  seen  for  300  kA. 

This  apparent  imderestimation  of  internal  stresses  should  not  be  entirely  unexpected,  as  the 
assumption  of  filamentary  conductors  is  certain  to  result  in  poor  estimations  of  the  field  strength 
in  regions  close  to  the  filaments  themselves.  In  particular,  the  effect  of  shielding  by  the 
conductor  sections  is  completely  omitted.  The  "pinching"  of  the  field  lines  caused  by  the 
exclusion  of  field  from  the  conductor  sections  may  well  increase  the  field  between  adjacent 
lengths  of  conductor  by  nearly  a  factor  of  2  in  this  design,  as  the  gaps  through  which  the  field 
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Figure  8.  The  Magnitude  of  the  Measured  Magnetic  Fidld  Compared  to  Simulations.  Note 
That  the  Complex  Pattern  Is  Largely  Due  to  the  Presence  of  the  Feed  Cable. 


5.8  tons  100  ksi 

copper-b  eryllium 
yield-strength  = 


♦  ♦  155  ksi 

Figure  9.  Calculated  Forces  on  One  Section  of  a  Square  Pancake  Coil. 
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must  thread  only  cover  50%  of  the  coil's  area.  Since  the  exact  extent  of  this  effect  is  not 
calculated,  internal  stress  estimations  can  only  be  described  as  semi-quantitative.  However,  this 
information  is  still  quite  useful  for  engineering  purposes. 

3.  Conclusion 

The  three  models  described  here  have  proved  vital  in  the  development  of  an  understanding  of 
the  key  issues  in  induction  launcher  design.  Based  on  the  results  of  modeling,  certain  rules  of 
thumb  have  been  developed: 

(1)  For  the  best  coil/plate  coupling,  the  plate's  area  should  be  equal  to  or  greater  than  the 
area  of  the  magnet's  bore.  This  guarantees  that  most  of  the  field  lines  will  be  forced 
to  bend  behind  the  plate,  where  they  do  the  most  good.  Plates  smaller  than  the  bore 
will  have  field  lines  bending  around  their  sides,  so  that  there  are  forces  "squeezing" 
the  two  sides  of  the  plate  together.  As  these  forces  are  at  right  angles  to  the  direction 
of  travel  and  tend  to  cancel  in  any  event,  they  represent  wasted  effort.  The  use  of 
small  plates  explains,  in  large  part,  the  poor  efficiency  of  the  box  coil  and  round 
pancake  coil  designs.  The  square  pancake  coil  design  uses  plates  with  the  same  area 
as  the  coils.  One  price  to  be  paid  for  this  improvement  is  a  somewhat  weaker 
structure,  as  the  slot  that  holds  the  plate  offers  less  support  against  forces  attempting 
to  squeeze  it  shut. 

(2)  Coil  conductor  sections  closer  to  the  plate  yield  the  highest  fields  and  forces,  so  it  is 
best  to  have  as  many  near  the  plate  as  possible.  This  suggests  that  a  pancake  coil 
design,  where  all  sections  are  equally  close  to  the  plate,  might  be  superior.  A 
serpentine  coil  also  offers  this  level  of  proximity,  but  generates  relatively  small 
fields. 
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This  latter  comment  raises  some  questions  concerning  the  relative  performance  of  the  box 
coil,  round  pancake  coil,  and  square  pancake  coil  designs  used  in  these  studies.  The  nearly 
identical  performance  of  the  box  coil  and  round  pancake  coil  can  be  explained  by  the  fact  that  the 
box  coil  had  a  conductor  section  running  directly  behind  (and  actually  in  contact  with)  the  plate. 
The  proximity  to  the  rear  of  the  plate,  along  which  the  largest  eddy  currents  flow,  explains  why 
the  box  coil  performed  so  well;  sections  further  from  the  plate  added  much  less  to  the 
performance.  In  addition,  the  pinching  effect  on  the  field  between  this  section  of  conductor  and 
the  plate  may  amplify  the  launch  forces.  In  the  pancake  coil  designs  described  here,  there  were 
no  equivalent  conductor  sections,  as  the  bus  bars  connecting  the  upper  and  lower  coils  ran 
vertically  down  one  comer.  The  latest  pancake  coil  designs  now  have  a  bus  bar  that  mns  down 
one  comer  behind  (and  in  contact  with)  the  plate  and  then  down  to  the  bottom  coil.  This  is 
expected  to  improve  performance  considerably. 

Tbe  explanation  for  why  the  square  pancake  coil  design  outperforms  (in  terms  of  measured 
AT.)  the  round  pancake  coil  design  (even  when  bore-sized  plates  were  used)  may  be  more  subtle. 
One  explanation  is  based  on  the  idea  that  two  linear,  antiparallel  currents  elements  repel  one 
another  more  effectively  than  curved  elements  would  because  the  Lorentz  forces  along  their 
lengths  are  all  in  the  same  direction.  In  the  case  of  curved  elements,  some  components  of  the 
Lorentz  forces  would  tend  to  cancel  each  other  out,  wasting  effort.  Tlius,  it  may  be  the  fact  that 
the  square  pancake  coil  offers  more  linear  segments  that  are  parallel  to  the  major  eddy  current 
element  running  along  the  back  of  the  plate  that  explains  the  improved  coil/plate  coupling. 

In  any  event,  the  square  pancake  coil  design  would  seem  to  epitomize  these  ideas  concerning 
performance  improvement.  It  also  carries  the  added  benefit  of  being  naturally  compact,  so  that 
more  launchers  can  be  mounted  on  a  vehicle.  One  disadvantage  that  always  accompanies 
increased  performance,  however,  is  increased  intracoil  stresses.  Future  design  efforts  will  center 
about  containing  these  stresses  more  effectively. 
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